arXiv:1507.08147vl [astro-ph.GA] 29Jul2015 


Astronomy & Astrophysics manuscript no. ms ©ESO 2015 

July 30, 2015 


The near-infrared radiation background, gravitationai wave 
background and star formation rate of Pop III and Pop II during 

cosmic reionization 

Y. P. Yangi’2, p. Y Wang^’^ and Z. G. Dai^’^ 

* School of Astronomy and Space Science, Nanjing University, Nanjing 210093, China 

^ Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, Nanjing 210093, China 
e-mail: fayinwang@nju.edu.cn 

July 30, 2015 


ABSTRACT 

Context. The transition from Population 111 (Pop 111) to Population 11 (Pop II) stars plays an important role in the universe history. Due 
to a huge amount of ionizing photons generated by Pop III stars, they begin to ionize the intergalactic medium (IGM) at the early stage 
of reionization. Meanwhile, the feedback from reionization and metal enrichment changes the evolution of different populations. The 
near-infrared radiation background (NIRB) and the background of gravitational waves (SBGWs) from these early stars will provide 
important information about the transition form Pop III to Pop II stars. 

Aims. In this paper, we obtain the NIRB and SBGWs from the early stars, which are constrained by the observation of reionization 
and star formation rate. 

Methods. We study the transition from Pop 111 to Pop II stars via the star formation model of different population, which takes into 
account the reionization and the metal enrichment evolution. We calculate the two main metal pollution channels arising from the 
supernova-driven protogalactic outflows and “genetic channel”. We obtain the SFRs of Pop III and Pop II and their NIRB and SBGWs 
radiation. 

Results. We predict that the upper limit of metallicity in metal-enriched IGM (the galaxies whose polluted via “genetic channel”) 
reaches Zj-nt = 10 “^'^Zq at z ~ 13 (z ~ 11), which is consistent with our star formation model. We constrain on the SFR of Pop III stars 
from the observation of reionization. The peak intensity of NIRB is about 0.03 - 0.2 nWm“^sr“‘ at ~ 1/jm for z > 6. The prediction 
of NIRB signal is consistent with the metallicity evolution. We also obtain the gravitational wave background from the black holes 
formed by these early stars. The predicted gravitational wave background has a peak amplitude of 8 x 10“'’ at v = 158 Hz for 

Pop II star remnants. However, the background generated by Pop III.2 stars is much less than Pop II stars, with a peak amplitude of 
Dgiv - 1.2 X 10“‘' at V = 28 Hz. The background of Pop lll.l shifts to lower frequencies, and the amplitude of Ggh, for Pop lll.l 
stars shows a minimum value at v =: 10 Hz, due to the lack of gravitational wave signals from the stars with 140 Mq < M, < 260 Mq. 
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1. INTRODUCTION 


The first stars, known as Population III (Pop III), form at the 
end of the cosmic dark age and open the reionization era. Pop 
[II stars can have different characteristics compared to Pop I/II 
stars. Some authors have stated that Pop III can be very massive, 
on the order of 500 - 6QQMq (Abel et al. 2002; Omukai & Palla 
2001, 2003; Bromm & Loeb 2004), while others have claimed 
that the masses are lower, on the order of 30 - 6QMq or IOMq 
(H osokawa et al. 2011, 2012; Stacy et al. 2012). The initial con¬ 
ditions for Pop III star formation are entirely determined by 
the basic parameters of the cosmology (Tegmark et al. 1997; 
Yoshida et al. 2003; Bromm et al. 2009). According to the pop¬ 
ular “bottom-up” hierarchical structure formation in the stan¬ 
dard cold dark-matter (CDM) model, the hrst stars form in a 
sufficient amount of cold dense gas in a dark matter halo with 
mass > 10® ~ 10 ^Mq, which collapses at redshift z ~ 20 - 30 
(Bromm et al. 1999, 2002; Abel et al. 2002). At hrst, these stars 
contain no elements heavier than helium, viz., metal free, so- 
called Pop III. 1. The second generations of stars originate from 
the environment inhuenced by earlier star formation, dehned 


as Pop III.2, which exist in a pre-ionized region and are af¬ 
fected by previous generations of stellar radiation due to photo¬ 
ionization of ambient neutral hydrogen of Pop lll.l stars. Be¬ 
cause of supernova (SN) explosions and stellar winds, heavier 
elements are introduced into the intergalactic medium (IGM), 
gradually increasing its metallicity. On the other hand, some new 
galaxies inherit the metals from their low-mass progenitors. Af¬ 
ter the metal elements enriching to a certain threshold Zcut, the 
Pop I/II stars form (Bromm et al. 2001; Bromm & Loeb 2003; 
Schneider et al. 2002, 2006a; Mackey et al. 2003). 

There has not been any direct observation of Pop III stars. 
Fortunately, there are some observations to constrain the early 
stellar population and reionization, e.g. the near-infrared ra¬ 
diation background (NIRB), Gunn-Peterson trough (Fan et al. 
2006) and the background of gravitational waves (Marassi et al 
2009; Pereira & Miranda 2010). NIRB contains some impor¬ 
tant information of early universe, which arises from accu¬ 
mulated emission from the early galaxy populations with a 
large range of redshifts. Because photons lose energy in cos¬ 
mic expansion, the optical and UV radiation from the early 
stars during reionization should leave a signature in the ex- 
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tragalactic background light at NIR bands with wavelength 
less than a few micrometers. Although most high-redshift 
galaxies are below the limiting magnitude of current detec¬ 
tors, we could observe their redshifted cumulative spectrum, 
which would be present in any background emission in the 
NIR band (Barkana & Loeb 2000; Salvaterra & Ferrara 2006; 
Wyithe & Loeb 2006; Kistler et al. 2009; Bouwens et al. 2010; 
Robertson et al. 2010; Fernandez & Shull 2011; Munoz & Loeb 
201 1). On the other hand, since Pop II galaxies at z ~ 6 are al¬ 
ready observed with Hubble Space Telescope (HST) and will 
be in a short time studied in great details with James Webb 
Space Telescope (JWST), by studying the NIRB intensity, we 
can find many properties of the early stars, such as the star 
formation rate of high-redshift stars whose photons are red- 
shifted, the clustering properties (by studying the anisotropy 
power spectrum of the NIRB), the classification of Pop III stars 
(by studying the number of the Lyman alpha bumps), the tran¬ 
sition from Pop III to Pop II stars (by studying the evolution 
of Lyman bumps), the reionization and metal enrichment (as 
we discussed in this paper). However, determining the contri¬ 
bution of these early stars at high redshifts from observations is 
a hard task. There are two reasons; first, low redshift and very 
faint galaxies make a significant contribution to NIRB, which 
is not measured accurately. Second, the extragalactic NIRB is 
hard to distinguish from the brighter zodiacal foregrounds in 
the local matter, (e.g. interplanetary dust within the solar sys¬ 
tem), and the stars and the interstellar medium (ISM) of our 
Galaxy (Hauser & Dwek 2001; Kashlinsky 2005), which is two 
to three orders of magnitude larger than 10 nWm^^sr"'. Recent 
works suggest that the intensity of NIRB cannot be larger than a 
few tenths nWm^^sr^' (Kneiske et al. 2004; Stecker et al. 2006; 
Franceschini et al. 2008; Gilmore et al. 2009, 2012; Finke et al. 
2010; Kneiske & Dole 2010; Inoue et al. 2013). Except for mean 
intensity, the anisotropy power spectrum of NIRB could pro¬ 
vide direct information of these early stars (Cooray et al. 2004; 
Kashlinsky et al. 2004; Fernandez et al. 2010; Cooray etal. 
2012a). The Pop III epoch contains less projected volume than 
the ordinary galaxy populations (e.g. Pop II stars), leading to 
larger relative fluctuations, and their anisotropy power spectrum 
represents the clustering behavior of sources in the IR regime. 
However, for the ordinary galaxies (containing Pop II stars), their 
angular power spectrum is nearly a power law over a wide range 
of angular scales (Cooray et al. 2004, 2012a). 

The UV photons (A ~ 1000A) produced at z = 9 during 
reionization will be redshifted to the NIR band with wavelength 
A ~ 1/Lim. Thus the NIRB would depend on the evolution 
of the early stars during reionization. So far, many works 
have obtained some information about reionization, e.g. the 
optical depth for electron scattering, the end of reionization 
redshift, and the escape fraction of ionizing photons, which 
can be used to constrain the SFRs of the early stars(Yu et al. 
2012; Wang 2013), and further constrain the intensity of 
NIRB and SBGWs. Recently, there are many studies to cal¬ 
culate the contribution of stars and galaxies responsible for 
reionization and NIRB (Kashlinsky et al. 2002; Santos et al. 
2002; Salvaterra & Ferrara 2003; Cooray & Yoshida 2004; 
Cooray et al. 2004; Kashlinsky et al. 2004; Kashlinsky 2005; 
Kashlinsky et al. 2005; Madau & Silk 2005; Magliocchetti et al. 
2003; Fernandez & Komatsu 2006; Kashlinsky et al. 2007a; 
Thompson et al. 2007a,b; Fernandez et al. 2010; Cooray etal. 
2012a; Fernandez et al. 2012; Kashlinsky et al. 2012; 
Fernandz et al. 2013; Fernandez & Zaroubi 2013; Yue et al. 
2013). As a crucial part of the entire story of the cosmic history, 
reionization reflects the transformation of neutral hydrogen in 


the IGM into an ionized state, which is due to ionizing photons 
generated by the early stars and galaxies. In the process of the 
early star formation, the IGM is enriched with metals that are 
dispersed by the first SNe and stellar winds. This process might 
reflect the evolution of the stellar populations, although it is 
not understood when the transition from Pop III to Pop II takes 
place, due to an uncertainty in the “critical metallicity” Zc„( 
(Bromm et al. 2001; Jappsen et al. 2009a,b), and the redshift 
evolution of the IGM metallicity (Barlow & Tytler 1998; 
Ferrara et al. 2000; Greif & Bromm 2006) and some other metal 
pollution channels, e.g. “genetic channel” (Schneider et al. 
2006b; Trent! & Stiavelli 2009). With more Pop II star forming 
(although fewer ionizing photons per stellar), there are more 
ionizing photons to ionize the IGM, and the ionized bubbles 
gradually overlap, allowing the mean free path of ionizing 
photons to increase rapidly. The average volume fraction of 
ionized hydrogen in IGM increase rapidly. At last the entire 
universe is almost completely ionized (Barkana (fe Loeb 2001; 
Robertson et al. 2010; Bromm 2013). 

In this paper. We consider two main channels of met¬ 
als pollution: the supernova-driven protogalactic outflows 
(Ferrara et al. 2000; Furlanetto & Loeb 2005) and the “genetic 
channel” (Schneider et al. 2006b; Trent! & Stiavelli 2009). The 
former assumed that metals are enriched via supperbubbles re¬ 
sulting from supernova explosion in protogalaxies, and the lat¬ 
ter suggested that the new galaxies inherit metals from the 
lower mass progenitor galaxies. Moreover, we constrain the SFR 
through the observations and reionization, e.g. optical depth 
for electron scattering as measured by WMAP and Planck, 
and the redshift of the end of reionization with a possible 
range from 5 to 10. During reionization, the metal elements 
enrich in IGM. Pop III.l stars originate from a freshly col¬ 
lapsed halo, and Pop III.2 stars exist in a pre-ionized metal- 
free region and they form via hydrogen deuteride cooling. Pop 

II stars originate from the dark matter halos that have been 
polluted by metal enrichment. Thus the SFRs of these early 
stars would depend on the the hydrogen reionization fraction 
and pristine fraction (Greif & Bromm 2006; Wang et al. 2009; 
de Souza et al. 201 1). Following Fernandez & Komatsu (2006), 
this NIRB consists of several contributions: the continuum emis¬ 
sion from stars themselves, the series of recombination lines, 
the free-free and free-bound continuum emission from ion¬ 
ized gas or nebula, and the two-photon emission. Based on 
the constraint on the SFRs of the early stars, we can obtain 
the contribution of the first stars to the high-redshift NIRB 
during reionization, which is lower than the current observa¬ 
tion (the total luminosity 1 ~ 10 nWm^^sr ') (Santos et al. 
2002; Salvaterra & Ferrara 2003; Cooray & Yoshida 2004; 
Kneiske et al. 2004; Kashlinsky et al. 2005; Stecker et al. 2006; 
Kashlinsky et al. 2007a; Franceschini et al. 2008; Gilmore et al. 
2009; Finke et al. 2010; Kneiske & Dole 2010; Gilmore etal. 
2012; Kashlinsky et al. 2012; Inoue et al. 2013), due to the 
foreground pollution. On the other hand, as pointed out by 
Fernandez & Zaroubi (2013), the Lyman a emission from Pop 

III stars at high redshift could result in “bump” in spectrum of 
NIRB. This shape of Lyman a “bump” is determined by the tran¬ 
sition from Pop III to Pop II stars. 

In addition, these early stars are predicted to collapse into 
black holes (except in the mass range 140-260Mo, which die as 
pair-instability supernova (PISN) (Heger & Woosley 2002)), and 
expected to be the sources of the stochastic background of gravi¬ 
tational waves (SBGWs) (Schneider et al. 2000; Buonanno et al. 
2005; Sandick et al. 2006; Suwa et al. 2007; Marassi et al 2009; 
Pereira & Miranda 2010). Based on the SFRs, we calculate the 
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SBGWs that were produced by these early stars, including Pop 
II, Pop III.2 and Pop III.l stars. At present, some gravita¬ 
tional wave interferometers are operating in the frequency of 
10 - 3000 Hz, e.g. VIRGO, Laser Interferometer Gravitational- 
Wave Observatory (LIGO). In the future, next-generation grav¬ 
itational wave detection will open a lower frequency window, 
e.g. Laser Interferometer Space Antenna (LISA) covering the 
frequency range 10 - 0.1 Hz, and Big Bang Observer (BBO) 

operating in the range 0.01 - 10 Hz. The signals of gravitational 
waves will open a new window for the study of the cosmic tran¬ 
sition from Pop III to Pop II stars. 

The structure of this paper is organized as follows. In Section 
2, we outline the stellar models, including the properties of dis¬ 
tinct populations and their SFRs during the era of reionization. 
In Section 3. we give our reionization model by considering the 
transition from Pop III to Pop II. In Section 4, we calculate the 
spectrum of NIRB, which is contributed by the early stars before 
the end of reionization. In Section 5, the gravitational wave spec¬ 
tra produced by Pop III and Pop II stars are presented. In Section 
6, discussions and conclusions are given. In the whole paper, we 
assume a flat ACDM model with Q,„ = 0.27, Qa = 0.73, and 
Hq-1\ km s“' Mpc^'. 


2. THE STELLAR MODELS 


with the mass range from 5 to 150 Mq. As pointed out by 
Johnson & Bromm (2006), we use an intermediate stellar IMF 
between Pop III.l and Pop II for Pop III.2 stars. We adopt the 
IMF given by Larson (1998) with Mf = 150 Mq, and the mass 
range is from 5 to 250 Mq. So the mean mass of Pop IIL2 is 
~ 40Mo, which is consistent with the typical mass given by 
Yoshidaetal. (2007) and Hosokawa et al. (2011). The normal¬ 
ization is given by 



(iM./(M.) = 1. 


(3) 


The mean stellar mass of one population is 


M. 





(4) 


Here we use the results from Lejeune & Schaerer (2001) and 
Schaerer (2002) to calculate the main stellar parameters, such as 
the intrinsic bolometric luminosity LJ°^(M,), the effective tem¬ 
perature the main-sequence lifetime and the 

time-averaged hydrogen photoionization rate RniiMt). At first 
we define x = log[o(M,/Mo). 

For Pop III stars (Pop III. 1 and Pop III.2), the parameters are 
given by 


In this section, we consider the emissions from early stars during 
reionization. Following Fernandez & Komatsu (2006), we cal¬ 
culate the emission from two stellar populations. First, Pop III 
stars can be divided into the first generation stars (Pop III.l), 
which their formation only depends on initial conditions of the 
early universe, and the second generation stars (Pop III.2), which 
exist in a pre-ionized region and are affected by previous gen¬ 
erations of stellar radiation due to photo-ionization of ambient 
neutral hydrogen of Pop III.l stars. They form via hydrogen 
deuteride cooling (Johnson & Bromm 2006) with less masses 
and extremely poor metallicity Z < IO^^'^Zq. In the hydrogen 
deuteride cooling process, the gas temperature would reach that 
of the cosmic microwave background (CMB) within a Hubble 
time. When the gas temperature reaches the CMB limit, the 
Pop IIL2 stars with typical masses ~ 40 Mq form, which are 
smaller than the typical mass of Pop III.l stars. For simplifica¬ 
tion, we assume that the Pop III.2 stars have a lower-mass dis¬ 
tribution, but have the approximate properties of the Pop III. 1 
stars (Ohkubo et al. 2009), including the intrinsic bolometric lu¬ 
minosity, the effective temperature, the main-sequence lifetime, 
and the time-averaged hydrogen photoionization rate. Second, 
the Pop II stars, as metal-poor stars with metallicity Z = 1 /50 Zq, 
form in clouds that can undergo metal and dust cooling. The 
clouds are able to fragment into smaller masses, leading to typi¬ 
cal masses less than those of Pop III stars. 


logio(LfVio) 

logio(T:-^^/K) 

logio(G/yr) 

logio(^/r//s"') 


0.4568-H 3.897.r-0.5297.x^ 

3.639 -H L501x - 0.5561x2 0.07005x2, 

9.785 - 3.759X -h 1.413x 2 - 0.186x2, 

f 39.29 -H 8.55x (SMq < M, < 9Mq) 

1 43.61 -t- 4.90x - 0.83x2 (9^^ < M. < 500Mo 


For Pop II stars, the parameters become 


logio(Lj“V-Lo) 

iogio(r:^'^/K) 

logio(T./yr) 

logio(^H//s^') 


0.138-h4.28x- 0.653x2, 

3.92 H-0.704X-0.138x2, 

9.59-2.79XH-0.63x2, 

27.80 -H 30.68X - 14.80x2 h- 2.50x2. 


(6) 


This stellar model gives the number of ionizing photons emitted 
per stellar baryon 


/ RHi{MZ)Teff{MZ) \ nip 

- W. -/—■ 

where nip is the proton mass, Y - 0.25 is the mass fraction of 
helium, and TeffiMY) is the effective stellar lifetime, given by 
Eq (51) in Section 4. The effective stellar lifetime Tejf may be 
less than the real lifetime t*, because some low-mass stars are 
still not dead. 


2.1. Pop III & Pop II 


2.2. The star formation rate 


We adopt the similar initial mass function (IMF) as same as 
Cooray et al. (2012a) for Pop III.l and Pop II stars. The IMF 
of Pop III.l stars is (Larson 1998) 


f(M,) oc m;' 



( 1 ) 


where = 250 Mq, and the mass range is from 5 to 500 Mq. 
For Pop II stars, we adopt the IMF given by Salpeter (1955) 


f(M,) (X m;2-35^ 


( 2 ) 


Following Greif & Bromm (2006), we consider that the Pop III. 1 
stars form in the minihalos with the minimum virial temperature 
T,,ir~ 102K.PopIII .2 stars form in metal-free halos above Tv/r = 
10"^K via HD cooling. When star formation begins, a sufficient 
amount of cold dense gas accumulating in a dark matter halo 
is needed. Because of the collapse of the dark matter halos, the 
baryonic gas becomes viral equilibrium with the dynamically- 
dominant dark matter halos. 


GMi, 

Rvir 


( 8 ) 
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where M;, is the dark matter halo mass, /?,,,> is the virial 
radius, and is the virial velocity. The virial radius is 
(Barkana & Loeb 2001; Bromm 2013) 


Ryir ~ 0 . 2 kpc 


M,, 


10^ Mp 


1/3 




(9) 


where A^. = PvirlPh is the overdensity after virialization nearly 
finishes, p,,,> is the virial density of dark matter halos, pt - ^hPcr 
is the baryon mass density of the universe, and ~ IStt^ in the 
Einstein-de Sitter model. The gas heats up resulting from the 
collapse of the dark matter halo, in which the virial temperature 
of the gas corresponds to the virial velocity of dark matter halos. 
So ksTpir ~ piUpV^,-!., which is leading to 


lO^Kp 


Ml, 


10^ Mp, 


2/3 



( 10 ) 


where p - 1 . 2 , 0.6 is the mean molecular weight for neutral and 
ionized primordial gas, respectively. For a dark matter halo with 
mass of IO^Mq, the virial temperature in this minihalo is ~ lO^K, 
which is below the threshold ~ 10 "^K of the cooling temperature 
of atomic hydrogen, so that the gas is unable to cool, resulting 
in no star formation. The gas would simply persist in hydrostatic 
equilibrium. However, the cooling in such a low-temperature pri¬ 
mordial gas could rely on molecular hydrogen (H 2 ) instead. The 
main formation channel is the sequence; H -t- e —> H + 7 , fol¬ 
lowed by H -H H —> H 2 -H e“ (Bromm 2013). At last, the evolu¬ 
tion of gas inside minihalos, driven by H 2 cooling, leads to the 
formation of Pop III.l stars with typical masses ~ lOOM©. After 
the Pop 111. 1 star formation, the composition of the primordial 
gas becomes different, and the cooling process would be more 
complicated. The hydrogen deuteride (HD) molecule provides 
an additional cooling channel. For such a dark matter halo with 
mass of > lO^Mo, the cooling of HD would be efficient with 
temperature ~ 10"^K, leading to the formation of Pop III.2 with 
the typical masses ~ TOM©. 

In order to obtain the SFRs of the early stars, we consider 
a semi-analytic approach by using the collapse fraction function 
of a dark matter halo. The famous one is the Press-Schechter 
formalism providing a way to calculate the abundance of the 
mass of dark matter halos (Press & Schechter 1974). For a given 
power spectrum P(k) oc \6k?', we adopt the matter power of 
Komatsu et al. (2009). The Gaussian variance of the fluctuations 
on the mass-scale M/, is 

Jo (11) 


where Mi, = {Aiy)np,,irR^j^ and W (k,Ri,) is a top-hat filter func¬ 
tion 


W = -^ [sin (kRi,) - (kRi,) cos (kRi,)]. (12) 

{kRhf 


The comoving number density of dark matter halos per unit mass 
could be given by Press & Schechter (1974) formalism 


_ ^mPcr /7 6c -5^/2a'[^^D^(z) ^IttCT^^ /I "23 

dMk Ml y 7T crMi,D (z)‘^ dlnMi, 

where the critical density of the universe is per = 1.8785 x 
10 “®/ 2 ^gcm“^, and 6c - 6p/p = 1.686 is the critical overdensity 
for a spherical perturbation. For the ACDM model, the growth 
factor is given by D{z) = g (z) / (1 +z)g (0), where 


with 


^m{z) 


Sl,„{\+zf 
flffl(l + z)^ + Qa 


(15) 


^^a(z) = 


_f^A_ 

flm(l + Z)^ + 


(16) 


A fraction of mass in the universe collapsing into halos with the 
mass more massive than Mk^nin, referred to as the collapse frac¬ 
tion, is written as 


2 pMi, 

Fcoi = — I Mi,nps(Mi„z)dMh, (17) 

Pm J 


where p„, - Sl,„pcr is the matter density of the universe, and 
nps {Mi„z) is obtained by the PS formalism. Mip^in and Mi,j„ax 
are determined by the virial temperature. Thus the preliminary 
SFR is given by (Greif & Bromm 2006) 


<A(2) = f*Pm^ 

iltn 


dF (;qI 

dz 

dz 

dt 


(18) 


where /, is the star formation efficiency. 

The transition from Pop III to Pop II stars is dependent 
on the evolution of metallicity. Here, we consider two chan¬ 
nel of metal pollution: supernova-driven protogalactic outflows 
(Ferrara et al. 2000; Furlanetto & Foeb 2005) and the “genetic 
channel” (Schneider et al. 2006b; Trent! & Stiavelli 2009). At 
first, we consider that metals are enriched via supernova-driven 
protogalactic outflows. For a protogalaxy within a dark matter 
halo, metals are blow out by supperbubbles (SBs) resulting from 
supernova explosion. According to Ferrara & Tolstoy (2000), the 
condition of blowout is that the blowout velocity v/, should be 
larger than the escape velocity Vg, which gives the fraction of 
the mechanical energy of superbubble that can blowout. The ef¬ 
ficiency of the metal escape produced by SNe is close to unity 
when blowout does take place (MacFow & Ferrara 1999). If a 
primordial galaxy has an exponentially stratified density distri¬ 
bution p oc exp(-z///), then the shock wave from SN explo¬ 
sion initially is decelerated by ISM and subsequently accelerated 
due to blowout in a lower density environment and continued 
SB luminosity. The velocity of shock wave has a minimum at 
z — 3H, which is defined as the blowout velocity. It is given by 
(Ferrara & Tolstoy 2000) 


Vb ^ 2.7 


pniptiH- 


1/3 


(19) 


where L is the mechanical luminosity of SBs, n is the num¬ 
ber density of a uniform ambient medium of protogalaxies, 
and H is taken to be a free parameter. For ISM parameters 
n - 0.5cm^, p - 1.25, blowout occurs if Vb > Vg. Thus, the 
critical mechanical luminosity should be 

IGM 

Lc^0.05pmpnH^\^-^\ . ( 20 ) 


The total mechanical luminosity L, must be larger than Lc when 
blowout occurs. Here we assume the escape velocity is equal to 
the circular velocity of the halo v'e ~ Vc - (GM;,/R/,)*^^. The total 
mechanical luminosity of SB in a galaxy could given by (Ferrara 
1998) 


giz) = 


_(5/2)a„(z)_ 

- Ga + (1 + a„(z)/2)(l + Qa(z)/70)’ 


(14) 


L, 


- SsNVftfb 


Q.b Mb 
fff (z) 


( 21 ) 
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where esn - 10^'erg is the characteristic energy of a SN, 
the cooling fraction of baryons fh - 1 in the halo 10"'< 
Tyir < 10^''K (Madau et al. 2001), and /, is the star forma¬ 
tion efficiency. f//(z) = (4;7rGp„, is the free fall time 

of dark matter halo. A typical mechanical luminosity is L, a; 
2.3 X lO-'^erg s“' for M;, = IO^Mq and z = 10. The mass in stars 
per SN event which is given by 

y = -■ ( 22 ) 


with Ml, ~ 10®Mq. Similarly, the effective mechanical luminos¬ 
ity leading to blowout is 

Lhiz, > Le) = Nma. (26) 

where Nc - LcIqbIssn is the number of SN in a cluster with me¬ 
chanical luminosity. Thus the fraction of the mechanical energy 
that could blowout is defined as (Ferrara et al. 2000) 

Lb max ILJc) 

Vb = — = —riTj -■ (27) 

Lt InNmax 


Here we assume the lowest mass of SN progenitors is 8Mq. Eq 
(21) gives the total mechanical luminosity. However, only a frac¬ 
tion of the mechanical luminosity could blowout from galax¬ 
ies, because the SNe might be occur in different regions in one 
galaxy, and some superbubbles with low mechanical luminosity 
could not blowout in certain regions. Thus we calculate the effi¬ 
ciency of blowout following Ferrara et al. (2000). Pop III.l stars 
are formed in minihalos, and recent simulation (Bromm et al. 
2009) implies that a minihalo only contains one Pop III.l star. 
Thus the star formation may be confined to a small region, lead¬ 
ing to SN explosions in the formation of a single superbubble. 
In this case, the single superbubble is likely to blowout L, > Lc 
and the fraction of the mechanical energy that could blowout is 
ijB - I- For larger galaxies where Pop III.2 and Pop II stars 
formed (Bromm 2013), the SNe would occur in different OB as¬ 
sociations, which are more widely distributed within galaxies. 
The luminosity function of OB association is approximated by 
(Oey & Clarke 1997; Portegies Zwart et al. 2010) 

^ ^ = KN-\ (1 < N < Nmax), (23) 


where N is the number of SN in a cluster. The probability to 
have N SNe in one OB association is p ^ and the average 
probability isp~ lINmax- The average number of SN per OB 
association Nsn ~ N,otlN,„ax, The total number of supernovae 
could be given by N,ot - ftfbi^hlLI,,,)Mi,lt // (z). The total num¬ 
ber of OB associations is K ~ Nmax- The mechanical luminosity 
of an OB association is Lqb - NssnIIob, where toB - 40 Myr is 
the time at which the lowest mass SN progenitor expires ~ SM©. 
The total mechanical luminosity could be given by 

^Ninax 

L,{z) = LoBdNoB = N,„ax ^ In Nmax - N,o, ^ - (24) 

Jl toB toB 

The first total mechanical luminosity given by Eq(21) is calcu¬ 
lated via the star formation rate of a proto-galaxy, which is re¬ 
lated to the mass of the dark matter halo. The second total me¬ 
chanical luminosity given by Eq(24) is to obtain the efficiency 
of blowout, which is related to the number of SNe in a cluster. 
Physically, the latter mechanical luminosity must be equal to the 
former. This relation between the maximum number of SN and 
the total number of SN is approximately 


Ntot — N„,ax In Nmax - 


(25) 


The above equation describes the dependence of the total num¬ 
ber of SN Ntot as a function of the maximum number of SN 
in a cluster, and we assume that the maximum possible num¬ 
ber of SN in a cluster is Nmax ^ 500, which is consistent with 
the Monte Carlo simulation of Ferrara et al. (2000). We find that 
Eq(25) could be applied to the case of Pop III. 1 stars in minihalos 


After the shock wave of SBs propagates to IGM, the metal 
bubble would be enriched to a larger zone due to the lower den¬ 
sity of IGM. The outflow will be confined by the pressure of 
IGM, which determines the radius of the metal bubble. The shell 
growth stalls when P, = Pq, where Pq is the IGM pressure in the 
surroundings of the galaxy as Pq = n{z)kBT, where n{z) is the 
average baryon number density of the cosmology, and the gas of 
IGM is heated to P ^ 2 x lO^'K by photoionization heating. Dur¬ 
ing the ionized zone overlapping, the SN-driven bubbles would 
propagate in the photoionization gas. Here we assume that the 
outflow propagates in ionized zone. 

We consider the standard evolution for an adiabatic, 
pressure-driven superbubble. The growth of the SN shell radius 
is (Weaver et al. 1977) 


254 pmpiiicM 


1/5 


(28) 


where njcM ~ n (z) is the number density of the ambient medium 
in IGM. The interior pressure is 




(38507r)2/5 




(29) 


When the shell growth stalls P,- = Pq, we have 
5 7^/^ 

..^ Ly^ipmpnbCMy^V"’ ( 30 ) 

(SSUtt) 


where - rjBLt is the effective mechanical luminosity, which is 
the fraction available for blowout into IGM. The stall age is 




73/4 

(550^ 


(pniptijCM) P\ 


3/4p-5/4,l/2 


0 


(31) 


The typical parameters are Pj » 33kpc and fj 1.5 Gyr for 
Mb - and z = 10. If z » 1 {LIm » Ga), the stall redshift 

is given by (Voit 1996) 





Hptst j 


(32) 


where G, = tism + ts- Lsmis the time during which a SN shell 
spreads through a galaxy, which is much less than the stall age in 
the IGM. The comoving metal-enrich radius is Pg ^ (1 -(- Zs)Rs- 
If C < I/Hq, the fraction of space with metals is approximately 
given by 



dMi,—Rl iMi„z,) 


drips 

dMb 


(Mh,z). 


(33) 


This equation could be well approximate if the metal-enriched 
bubble did not overlap and the expanding time is much less than 
the Hubble time //^*(z) . If the protogalaxies are randomly dis¬ 
tributed, then the filling factor would be p'^iz) = 1 - exp[-2'(z)]. 
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In fact, due to clustering, some new halos form in the metal- 
enriched regions. The fraction of space with metals should be 
corrected by an excess probability that two galaxies are located 
near each other, viz, the galaxy two-point correlation function 
^gg(Re) = bobm^hhiRe) (Furlanetto (& Loeb 2005), where is 
the correlation function of dark matter halo (e.g. Greif & Bromm 
(2006)), bQ ^ b{Mh,miii) is the bias of the newly formed galaxies 
and b,„ is the bias of the metal enriched regions, which is given 
by 

f dMh{An l3)pRlb(m)(dnps I dMh) 

bm =--■ (34) 

J dMhMhidnps IdMh) 

Thus, the probability that a new halo lies within metal-enriched 
region is then approximately given by 

Qe^Q'A^+^ggiRe)] (35) 


Finally, the filling factor would be Pe{z) = 1 - exp[-2e(z)] after 
assuming the wind hosts are randomly. 

On the other hand, the metal enrichment of galaxies can pro¬ 
ceed via “genetic channel” (Schneider et al. 2006b), that is, met¬ 
als are enriched via the merger of the lower mass progenitors, 
rather than through outflows from neighbours. Here, we calcu¬ 
late the probability fouiz) that a new collapsing halos accreting 
onto some old halos via Extended Press-Schechter model (the 
Appendix of Furlanetto & Loeb (2005)), which is 

nOO 

L dAli^AIfj(^dnps jdAIi2)F(^<: Aliifyiif^^Zh\^h->z) 

foMiz)^^ -, (36) 

dMhMhidnps!dMh)F{< Mh,mm,Zh\Mh,z) 

where the factor 2 is attributed to the assumption that a parent 
halo with a mass Mu < 2Mhj„i„ would be included in the new 
halo component, and the fraction of the accreted mass in a halo 
with a mass M;, at redshift z is 


F(<Mh ,min ■> Zh\Mh,z) = erf 


Sc(zh) - 6c{z) 




Mhjnin ^Mk ' 


(37) 


where Zh corresponds to some earlier time when a parent halo 
with a mass (< formed, which is fixed by the dynamic 

time within a galaxy (Furlanetto & Loeb 2005). After a halo 
merging a mass above it would most likely form Pop II 

stars. Thus, the probability Pgdz) that a dark matter halo did not 
inherit any metals from its progenitor is Pgdz) - 1 - fouiz)- 

In Ligure 1, the solid curve corresponds to the evolution of 
the pristine fraction Ppns(z) = 1 - Peiz), and the dashed curve 
corresponds to the evolution of Pgdz)- We And that the metals of 
IGM begin to be significantly enriched via supernova-driven pro- 
togalactic outflows at z ~ 10. However, due to Pgdz) ~ 0.2 - 0.6 
during a large range redshift, the metal enrichment of galaxies is 
still dominated by “genetic channel” at z > 6. 

Since the Pop III stars generate many ionizing photons, the 
ionized bubbles would be photoheated to ~ 10"^K, which prevent 
the ionized gas to collapse fresh stars. Thus, the formation of Pop 
III. 1 stars is suppressed by a factor being equal to the volume 
Ailing fraction of ionized regions Qiodz), which is discussed in 
the next section. The SLR of Pop III. 1 stars is given by 


^m.dz) 


X 


/.,///.lPm^(l -Pdz))pgdz) 

i 4»»j 


(1 - QUz)) 


coll 

7',„v=10‘'K 

dz 

dz 


dt 


(38) 


As pointed out by Johnson & Bromm (2006), the free electrons 
can boost the production of H 2 , leading to a lower tempera¬ 
ture where the HD can be cooled. These stars, originating from 
metal-free gas and cooling via HD channel, would be less mas¬ 
sive than Pop III. 1 stars. The formation of Pop III.2 stars requires 
an increased abundance of free electrons. Here we consider two 
main pathways toward Pop III.2 (Bromm 2013). The first one 
results from the photo-ionization of ambient neutral hydrogen of 
Pop III. 1 stars, and the non-equilibrium recombination leads to 
a boosted abundance of H 2 and HD after the Pop III stars have 
died (Yoshida et al. 2007). In order to prevent pre-enrichment 
of the gas, these Pop III stars (260 < m < 500Mo) have to 
directly collapse to black holes, the number fraction of which 
is f„ - dMtf(Mt) ^ 6 x 10“^. The second one arises 

from the collision ionization in shocks that originate from the 
collapse of metal free gas of more massive dark matter halos 
with = lO^K (Greif & Bromm 2006). Thus, the SLR of Pop 
III.2 is 


Q/, 

^pni.liz) = - Pdz))Pgdz){f*,IIIAfnQwn{z) 

Si Im 


dF(;Qll 


dz 


= 10'‘K 


= 103K 


+ /*,7/7.2 


dFcoii 


dz 


= 10-‘K 


■(39) 


In fact, the first term in the above equation is small due to the 
number of the stars that directly collapse to black hole is very 
few. Lor Pop II stars, we assume that they are formed from the 
metal-enriched dark matter halos with Ty,, = 10"^K. The SLR of 
Pop II is 


^idz) = /.,77Pm^[l - (1 -Pdz))pgdz)] 

IZii) 


dF coll 



dz 

dz 

Tvi 

= 10‘‘K 

dt 


. (40) 


In Eq(38) and Eq(39), the SLRs depend on the hydrogen reion¬ 
ization fraction Qion that is as a function of redshift z, which is 
discussed in the next section. 

Due to inhomogeneous metal pollution, some regions where 
the metal enrichment is overlapped have larger metallicity. The 
upper limit of metal-enriched IGM metallicity in these rich re¬ 
gions could be estimated by 


^rich 


(z) = 


yn 


Peiz)pb 


r 


rjh,min,z'^^e(.Z ^4^11 (z ) 


dt 

dz' 


dz', (41) 


where the metal yields are yu - 0.005 for Pop II, which is con¬ 
sistent with the proposed values in Greif & Bromm (2006), pb 
is the mean baryon mass density of the universe, ^dz)4^n{z) is 
the SLR of Pop II stars which are polluted via supernova-driven 
protogalactic outflows, and ^dz) - Pe{z)l[l - (1 - Pe{z)pgdz))]- 
The metals in the overlapped metal-enriched region are mainly 
from the Pop II stars which are polluted via SNe, thus we ignore 
the first contribution of metal enrichment via the SNe of Pop III 
in this overlap region. On the other hand, we calculate the up¬ 
per limit of the metallicity of the galaxies whose all progenitors 
were polluted by “genetic channel”, which is 


Zgc (z) = 


yn 

PiSM 


r 


7fy(l - rjB(Mh,mm,z')Xgdz')i/^n(z') 


dt 

dz' 


dz', (42) 


where (gdz)4^ii{z) is the SLR of Pop II stars which are pol¬ 
luted via “genetic channel”, and f,gdz) = [1 - Pgdz)]l[^ - (1 - 
Pe(z)Pgdz))]- We assume that the mean baryon mass density of 
ISM is pisM - 10“^^cm“^, and the ratio of the mean distance 
between two galaxies to the scale of a galaxy is tc,, ^ 100. 
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3. THE REIONIZATION 


The volume filling fraction of hydrogen ionized regions Q'.^^ is 
(Madau et al. 1998) 


^Q'ion _ V /. . _ 9^ 

dt Hh.Q ^ tree ’ 


(43) 


where i represents Pop II, Pop III.2 and Pop III.l, fesc is the 
escape fraction of ionizing photons. In fact, due to clustering, 
the probability that a fresh halo lies within an hydrogen ionized 
(shown in Eq(38) and Eq(39)) is 


Qion - 2;on[i + ^hh{z)] 


(44) 


Irec is the volume averaged recombination time, which is given 
by 

'tree = [CH//(zX‘^nH.o(l + z)\l + T/4Y)]-\ (45) 

whereCH//(z) = is the clumping factor of ion¬ 

ized hydrogen. We use a simple analytic fit of the form Chii (z) - 
1 + 9[(1 +z)/7]-^ for z > 6, and Ch//(z) = 10 for z < 6 
(Greif & Bromm 2006). We assume that the mass fractions of 
hydrogen and helium are Y = 0.75 and V - 0.25, respectively. 
Thus the mean hydrogen number density at z = 0 is given by 

riHo- ~ 1.956 X 10^’cm"^. (46) 

nip 


The number of hydrogen-ionizing photons per stellar mass dur¬ 
ing stellar lifetime is 


qi = 


Ml 


(47) 


where i presents Pop II, Pop III.2 and Pop III.l stars, and 
TeffiMt) is the effective stellar lifetime, which may be less than 
T,, given by Eq. (51). The optical depth for Thomson scattering 
determined by the ionization history is 


cctt f dz'Q,on(z')nff,o(l+z'f(l+Y/4X)-^. (48) 

Jo dz' 


In order to constrain our star formation model with fewer 
parameters, we adopt a fixed optical depth of Thomson scattering 
T - 0.08, which is satisfied with the observation of WMAP nine- 
year data with t = 0.089 + 0.014 (Hinshaw et al. 2013) and the 
observation of Planck results with t = 0.066 + 0.016 (Ade et al. 
2015). The star formation efficiency of Pop II stars is assumed 
to be /».// = 0.01, which is constrained by the observation of 
the SER at z ~ 5 - 10 (Bouwens et al. 2012a,b; Schenker et al. 
2013). We take the same value of the star formation efficiency 
ft.iii for both Pop III.l and Pop III.2 and the same value of the 
escape fraction of ionizing photons for all populations fesc- We 
assume that stars form from redshift Zm - 30. Eigure 2 shows 
the three epochs of the cosmic star formation history: the blue 
line denotes Pop II stars, the black line denotes Pop III. 1 stars, 
and the red line denotes Pop III.2 stars. Because the negative 
feedback from a star forming in a dark matter halo prevents the 
formation of other stars in the same halos, leading to Pop III stars 
have a lower star formation efficiency than Pop II stars, and we 
set ft III < ft II. Eor the same optical depth t = 0.08, we consider 
three cases: Case A (/,/// = O.Ol, fesc - 0.21), Case B (ftiii - 
0.004, fesc = 0.45) and Case C (/.,/// = 0.002, fesc = 6.75), 
which are denoted by dashed, solid and dotted lines, respectively. 


Eor the same t value, a larger value of fesc will result in a lower 
SER and a lower value of fesc in a larger SER. We also find that 
the SER of Pop II raising earlier due to the metal pollution via 
“genetic channel”. 

Eigure 3 shows the reionization history for three star forma¬ 
tion models. The end of reionization is Zend - 5.1, 6.9, 7.2, cor¬ 
responding to Case A, B and C, which are denoted by dashed, 
solid and dotted lines, respectively. We find that larger SER of 
Pop III. 1 stars causes a larger hydrogen reionization fraction at 
high redshift, however, lower fesc leads to a later end of reion¬ 
ization. Due to larger SER of Pop III. 1 at high redshift and a 
huge amount of ionizing photons generated by Pop III. 1 stars, 
the hydrogen reionization fraction Qion raises at a higher redshift 
z ~ 20 until the SER of Pop III. 1 becoming low at z ~ 10. Soon 
the Pop III.2 and Pop II stars rapidly form, leading to the hydro¬ 
gen reionization fraction Qion raising rapidly until the universe 
completely ionized. 

According to Eq (41) and Eq (42), we obtain the evolution of 
the upper limit of the metallicities in the metal-enriched region 
of IGM and in the galaxies whose all progenitors were polluted 
via “genetic channel”, as shown in Eigure 4. The upper limit 
is attributed to the inhomogeneous metal pollution. We predict 
that: (i) for the metal-enriched region of IGM, the upper limit 
of the metallicity reaches Z^,-, = 10^^ ^Zq (Bromm et al. 2001; 
Schneider et al. 2002) at z ~ 13; (ii) for the galaxies polluted via 
“genetic channel”, it reaches the critical value at z ~ 11. Pop 
II stars with low mass formed for Z > Zerit, which is in good 
agreement with the star formation model in Eigure 2. Note that 
as shown in Eigure 1, the metal enrichment is dominated by “ge¬ 
netic channel”, however, at a certain redshift, the upper limit of 
the metallicity of the metal-enriched region is larger than that of 
the galaxies polluted via “genetic channel”. The reason is that 
the fraction of space with metals is so small that the metallic¬ 
ity would be large, if the Pop II stars that were polluted via SN 
outflows always form in this metal-enriched region. 


4. THE NEAR-INFRARED BACKGROUND 


The intensity of the NIRB offers a window of probing the era of 
reionization. Eollowing Eernandez & Komatsu (2006), we cal¬ 
culate the NIRB from the epoch of reionization. There are sev¬ 
eral contributions to the emission of NIRB: the continuum emis- 
sion from stars themselves L^, the series of recombination lines 

, the free-free and free-bound continuum emission from the 

—cont — 2y 

ionized gas or nebula , and the two-photon emission Ly . 
Eor the first stars, one of the remarkable properties is a “bump” 
in the spectrum of NIRB, which is from the Lyman a emis¬ 
sion (Eernandez & Zaroubi 2013). The Lyman a bump would be 
higher if Pop III stars were more massive and presented at lower 
redshifts, and it would evolve with the transition from Pop III 
stars to Pop II stars. Eor the IGM, the hydrogen density is lower 
than that of the stellar nebulae. Here we neglect the emission 
from IGM, because it is only a small part of NIRB (Cooray et al. 
2012a). The intensity of the NIRB is given by 


j _ f , p((l +z)v,z) 

47r J ^ //(z)(l+z) 


(49) 


where the volume emissivity is given by 


p(v, z) = 

Mt 



dMtf(Mt)L,iMt)Teff(Mt). 


(50) 


The stellar effective lifetime could be approximately given by 
TeffiMt) = min[T.(M.), Ts (z)]. (51) 
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Ly(Mt) is the time-averaged luminosity in the frequency inter¬ 
val dv for a radiative process a for stellar or nebular component, 
which consists of the stellar blackbody emission, as well as the 
reprocessed nebular emission, such as two-photon emission, re¬ 
combination line, free-free and free-bound continuum emission. 
Ts is the time from the formation of the first stars of the universe 
to the age of the universe at redshift z (Fernandez & Komatsu 
2006). 

Following Fernandez & Komatsu (2006), we calculate four 
contributions to the emission of NIRB: (1). The stellar spectrum, 
which is the Planck function with the Lyman absorption: 




AnRl{M,)By{rJf{M,)), hv <13.6 eV, 
0, /2v>13.6eV, 


(52) 


where By(Teff^ is the Planck spectrum. R, is the stellar ra¬ 
dius which is determined by the intrinsic bolometric luminosity 
(M,) and the effective temperature (M.) are given by 
Eq(5) andEq(6), respectively; (2). The luminosity of two-photon 
emission, which is given by 


= —(1 - fLya)P{vlVLya)RHl{M,), (53) 

VLya 


where viya - 2465THz, fiya - 0.64, P{y)dy is normalized 
probability of generating one photon via two-photon decay in the 
range dy = dvjvLya for y = v/v^^q, < 1 (Eemandez & Komatsu 
2006); (3). The line luminosity, which is given by 

Ly\M^) = flyahVLya^iv - VLya)RHl{MP), (54) 


where (p{v - viya) is the line profile, which is taken to be a S- 
function (p(v-vrya) = 6^{v—viya). Note that the above equation 
is for the Lyman alpha line only. (4). The free-free and free- 
bound continuum luminosity, that is 


X 

where the case B recombination coefficient is ^ 2.17 x 
10-iOj'-o.7395, ^^^and g/bin) are the Gaunt factors for free-free 
and free-bound emission respectively, approximately, 'gff ^ 1.1 
and gfbin) = 1.05. Tg is the gas temperature, here we assume 
that Tg ss lO^K, and the line profile is 02(1O"^K) ^ 1.0. Then 
Rhi{M,) is given by Eq(5) and Eq(6). 

Eigure 5 shows the spectra of NIRB at z > 6. Case A, B 
and C are denoted by dashed, solid, dotted lines, respectively. 
Different color lines denote different contributions to the emis¬ 
sion of NIRB. The Lyman a and two-photons emission are dom¬ 
inant. We find that larger value of /,,/// leads to larger Lyman 
a bump and stronger two-photon emission. The mean intensity 
of NIRB contributed by these early stars during reionization is 
nearly ~ 0.03 - 0.2 nWm^^sr * at z > 6. We also plot the con¬ 
tributions of the different stellar populations to NIRB, as shown 
in Eigure 6. The blue, black and red lines denote the NIRB con¬ 
tribution of Pop II, Pop III. I and Pop III.2, respectively. We find 
that the component of the Pop III.2 dominates the total NIRB 
spectra and the two Lyman a bumps (~ 1215(1 -i-z)A) associated 
to Pop IIL2 and Pop III. 1 are obvious at z > 6. If we can de¬ 
tect the Lyman a bumps, it would reveal many properties of the 
early stars, such as the star formation rates, redshift distribution 
and the classification of the Pop III stars. 


T,—1/2 

, J- a — 


6.8 X 


a'' 


*//+15^ L —ns—ww 

* n=2 


(55) 


5. STOCHASTIC BACKGROUND OF 
GRAVITATIONAL WAVES 


In this section, we calculate the SBGWs that are generated by 
these stars collapse to black holes. The SBGW is mainly depen¬ 
dent on the SER and IME of Pop III and Pop II stars. The flux 
received in gravitational waves is 


Rviyobs) — 


f 


1 dEcw dv 

— - ; -tA(z) — dM^dV. 

dv dvobs M* 


47Tdj^ 


(56) 


where di is the luminosity distance, dEcwldv is the specific en¬ 
ergy of the source, is the star formation rate, fiM,) is the 
IMP of one population, and M, is the mean stellar mass given 
by Eq(4). As pointed out by Carr (1980), the specific energy flux 
per frequency is 


fv{Vobs) — 


1 dEcw dv nc^ ^ 

And]^ dv dvobs 2G 


(57) 


where hsH is the dimensionless gravitational wave amplitude. 
The flux of gravitation waves would be 

nc^ 2 

EyiVobs^ — (58) 

ZLr 


The integral dimensionless gravitational wave amplitude pro¬ 
duced from all events that stars collapse to black holes is given 
by (Pereira & Miranda 2010) 


/j2 = J_ r hi ^(z)^^dM,dV. (59) 

Vobs J M. 

The comoving volume can be expressed as 

dV = 4;r4 j [q,„ (1 -H z)^ -H dz. (60) 


Por the case that a star collapse to a black hole, we as¬ 
sume that gravitational waves radiate with an efficiency ecw = 
AEcwInirC^, where AEcw is the total gravitational waves energy 
and nir is the mass of black hole. The efficiency ecw ^ 7 x 10“^, 
if the collapse is axisymmetric (Stark & Piran 1985). The char¬ 
acteristic amplitude is given by (Thorne 1987) 


'T-BH 




di ^-1 


‘^"''Mo^'lMpc 


( 61 ) 


and the observed gravitation waves frequency is 


^ L3x104Hz(^)(1+z)-‘. 

Mr 


(62) 


We consider that the black holes are formed with M, > 25Mq 
for Pop II and Pop III. Por the stars with 25 < M, < 140Mo, 
the black holes have the same mass of the helium core of their 
progenitors (Heger & Woosley 2002) 


13 

niy = niHe = “ 20Mo). 


(63) 


In the range 140 < M, < 260Mo, the stars are completely dis¬ 
rupted in PISN explosions, leading to no black holes, niy - 0. Por 
M, > 260Mo, these stars would directly collapse to black holes, 
thus we neglect stellar mass loss and the masses of the black 
holes are equal to those of the their progenitor stars, niy - M,. 
According to Eq (62), the range of frequency is determined by 
the mass range of IMP after assuming that the largest redshift 


Article number, page 8 of 12 












Yang, Wang & Dai: NIRB, SBGWs and SFR of Pop III and Pop II during cosmic reionization 


of star formation is z = 30. For Pop II, the frequency is Vghs > 
6.4Hz; for Pop III.l, the background is shifts to lower frequen¬ 
cies due to the direct collapse for M* > 260Mo, Vghs > 0.6 Hz; 
and for Pop III.2, Vohs > 6.4Hz, which has the same range as 
Pop II stars, because the stars with 140 < M, < 25QMq were 
disrupted by PISN. 

The gravitational energy density parameter Slew is defined 
as the closure energy density per logarithmic frequency span 
(Pereira & Miranda 2010) 


Slew 


1 dpew 
Per d In 


3hI 


^Ibsh 


2 

Be¬ 


rn 


In Figure 7, Pop II, Pop III.2 and Pop III.l stars are denoted 
by blue, red, black lines, respectively. The sensitivity curves of 
advanced LIGO HlLl, LISA and BBO are denoted by green, or¬ 
ange and purple, respectively (Thrane & Romano 2013), assum¬ 
ing T = 1 yr of the observation. Case A, B and C are denoted by 
dashed, solid, dotted lines, respectively. The predicted gravita¬ 
tional wave background has a peak amplitude of Slew - 8 x 10“^ 
at V = 158 Hz for Pop II star remnants. However, the background 
generated by Pop III.2 stars is much weaker than Pop II stars, 
with a peak amplitude of Slew — 1-2 x 10“^* at v = 28Hz. The 
background of Pop III.l shifted to lower frequencies, and the 
amplitude of Slew for Pop III. 1 stars shows a minimum value at 
Vohs - lOHz, due to the lack of gravitational wave signals from 
the stars with 140 < M, < 260Mq. As shown in Figure 7, it is 
difficult to observe the SBGWs from the early stars for the gravi¬ 
tational wave detectors. However, as shown in Case A, we might 
have an opportunity to detect the SBGWs signal from Pop III.l 
stars at observed frequency Vobs - 3.5 Hz, which are twice than 
the sensitivity of BBO detector. 


6. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have constructed the star formation history for 
Pop III and Pop II via the collapse function of dark matter ha¬ 
los. These SFRs, calculated by the collapse fraction function, are 
self-consistent with the observations of reionization (the reion¬ 
ization optical depth measured by WMAP and Planck, reion¬ 
ization redshift in range from 5 to 10) and NIRB (the total lu¬ 
minosity 1 ~ lOnWm^^sr '). At first. Pop III stars form from 
pristine baryonic gas in dark matter halos. Due to the metal en¬ 
richment raising in the universe, more and more Pop II stars 
formed in metal-enriched region. There are two main channels of 
the metal pollution: the supernova-driven protogalactic outflows 
(Ferrara et al. 2000; Furlanetto & Loeb 2005) and the“genetic 
channel” (Schneider et al. 2006b; Trent! & Stiavelli 2009). The 
former is assumed that the new dark matter halos are formed 
in the IGM that is enriched with the metals dispersed by the 
first SNe and stellar winds, and the latter suggests that the new 
galaxies inherit metals from the lower mass progenitor galax¬ 
ies. After the metals in the IGM enriching to a critical threshold, 
the Pop I/II stars form gradually. Therefore, the SFRs of these 
populations would depend on the hydrogen reionization fraction 
and pristine fraction, which are shown in Eq (38), Eq (39) and 
Eq (40). Our results show that the metal enrichment is dominated 
by “genetic channel” during a large redshift range, leading to the 
earlier raising of the SER of Pop II stars, which agrees with the 
result of Trent! & Stiavelli (2009). However, the upper limit of 
the metallicity of the metal-enriched region could be larger than 
that of the galaxies polluted via “genetic channel” due to the 
inhomogeneous metal pollution of the supernova-driven proto¬ 
galactic outflows. 


Eor the case of the supernova-driven protogalactic outflow, 
we make a more reasonable assumption that the SN winds would 
stall if their pressures are equal to the pressure of IGM, which 
is different from the previous papers (Eurlanetto & Loeb 2005; 
Greif & Bromm 2006) that calculated the process of the propa¬ 
gation and distribution of metals by assuming that the SN wind 
propagated for half of the age of the universe via Sedov solution 
(Eurlanetto & Loeb 2005), however, Sedov solution is not well 
approximated after the SN winds have been stalled by IGM. Our 
result shows that the IGM is enriched via SN outflows during 
z = 4~ 10. 

In this paper, we study the NIRB, the SERs of different pop¬ 
ulations, and the reionization history simultaneously, which is 
different from Eernandez & Komatsu (2006) and Cooray et al. 
(2012a). In previous works, e.g. Cooray etal. (2012a) and 
Eernandez & Zaroubi (2013), the transition of Pop III to Pop II 
stars is given by fp - (1/2){1 + erf[(z - ZtllcTp]), where fp is 
the fraction of Pop III stars, Zt is the transition redshift and cTp 
is the length of the transition. The above equation describing the 
transition of Pop III to Pop II might be too simple. In fact, the 
transfer from Pop III to Pop II stars is a very complex process. 
Here, we consider the some main properties evolution of three 
stellar populations, such as the metallicity, reionization and star 
formation rate. Einally, these properties would affect the charac¬ 
ter of the spectra of NIRB. 

As pointed out by Eernandez & Zaroubi (2013), the Lyo- 
bump could reveal information of the Pop III era. Our results 
show that the Pop III would make a main contribution to the 
NIRB at high redshift and the spectra of the NIRB might 
show two bumps due to the different distributions of Pop 
III.l and Pop III.2 stars. However, the predicted intensity 
(< 0.2 nWm^^sr ') of NIRB from high redshift (z > 6) 
stars is much lower than the contribution of the foreground 
(~ 10 nWm^^sr^^), leading to measuring the signal from 

high-redshift stars is very difficult. The contributions from 
low-redshift galaxies and the bright zodiacal foreground should 
be correctly subtracted. Many works have attempted to measure 
the excess of NIRB without the contributions of the low- 
redshift galaxies and other foregrounds (Dwek & Arendt 1998; 
Gorjian et al. 2000; Kashlinsky & Odenwald 2000; Totani et al. 
2001; Wright 2001; Kashlinsky et al. 2002; Magliocchetti et al. 
2003; Salvaterra & Eerrara 2003; Cooray & Yoshida 2004; 
Kashlinsky et al. 2004; Kashlinsky 2005; Matsumoto et al. 
2005; Kashlinsky et al. 2007b; Thompson et al. 2007a,b; 
Cooray et al. 2012b; Kashlinsky et al. 2012). On the other hand, 
Eernandez & Zaroubi (2013) suggested that we can detect the 
relative change in the intensity of NIRB, which results from a 
Lyman a bump as a function of wavelength. However, there are 
still many ways causing such a change, such as the evolution of 
the escape fraction or a rapidly changing star formation rate as a 
function of redshift. 

We also calculate the stochastic background of gravitational 
waves from the collapse of the early stars. However, theo¬ 
retically, there are several astrophysical sources contributing 
to the background of gravitational waves, including transient 
sources (e.g. compact binary coalescence, supernovae, gamma- 
ray bursts etc.), long-lasting transient sources (e.g. magnetars, 
long gamma-ray bursts etc.) and continuous sources (e.g. pul¬ 
sar). Thus it is necessary to distinguish between different sources 
from the spectrum of gravitational waves. In this paper, we And 
that it is difficult to observe the SBGWs signals from the early 
stars, however, for /,_/// 0.01 and fesc - 0.21, BBO operat¬ 

ing in the range 0.01 - lOHz might detect lower frequency signal 
from Pop III. 1 stars, which are twice than the sensitivity of BBO. 
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Fig. 1. The evolution of the non-metal fraction. The solid curve cor¬ 
responds to the evolution of the pristine fraction PpriAz) = 1 - pAz), 
and the dashed curve corresponds to the evolution of PgAz) that a dark 
matter halo did not inherit any metals from its progenitor. 



Fig. 2. The three epochs of the cosmic star formation history; the blue 
line denotes Pop II stars, the black line denotes Pop III. 1 stars, and the 
red line denotes Pop III.2 stars. Case A (/.,;// = 0.01, /„c = 0.21), Case 
B (/../;/ = 0.004, fesc = 0.45) and Case C (/.,/;/ = 0.002, f^sc = 0.75) are 
denoted by dashed, solid and dotted lines, respectively. The observed 
data is given by Bouwens et al. (2012a,b) and Schenker et al. (2013) 


Fig. 3. The hydrogen reionization fraction Qia„ as a function of redshift 
z. Case A (/,,/;/ = 0.01,/„c = 0.21), Case B = 0.004,= 0.45) 
and Case C (/.,//; = 0.002, fe^c = 0.75) are denoted by dashed, solid and 
dotted lines, respectively. 



Fig. 4. The evolution of the upper limit of the metallicity. The soild line 
denotes the upper limit of the metallicity in the metal-enriched region 
of IGM. The dashed line denotes the upper limit of the metallicity in 
the galaxies whose all progenitors were polluted via “genetic channel”. 
The dotted line denotes the critical metallicity Z^i, = 10“^'^ Zg. 


Article number, page 11 of 12 









A&A proofs: manuscript no. ms 




log(v„,./Hz) 


Fig. 5. The spectra of NIRB at redshift z > 6. Case A = 

0.01,/psc = 0.21), Case B = 0.004,= 0.45) and Case C 
(/•,/;/ = 0-002, fcsc = 0.75) are denoted by dashed, solid and dotted 
lines, respectively. Different color lines denote different contributions 
to the emission of NIRB. 



Fig. 7. The spectrum of the gravitational waves. Pop II, Pop III.2 and 
Pop III. I stars are denoted by blue, red, and black lines, respectively. 
Case B (/.,/// = 0.004, f^sc = 0.45) and Case C (/.,;// = 0.002, fg^c = 
0.75) are denoted by dashed, solid and dotted lines, respectively. The 
sensitivity curves of advanced LIGO HlLl, LISA and BBO are denoted 
by green, orange and purple, respectively (Thrane & Romano 2013), as¬ 
suming r = 1 yr of observation. 


Fig. 6. The contributions of the different stellar populations to NIRB. 
The blue, black and red lines denote the NIRB contribution of Pop II, 
Pop III.l and Pop III.2, respectively. Case B (/.,;// = 0.004, = 0.45) 

and Case C = 0.002, fe^c = 0.75) are denoted by dashed, solid and 
dotted lines, respectively. 
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